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TECHNICAL NOTE 3^99 


CALCULATION OF THE SUPERSONIC PRESSURE DISTRIBUTION 
ON A SINGLE-CURVED TAPERED WING IN REGIONS 
NOT INFLUENCED BY THE ROOT OR TIP 
By Walter G- Vincenti and Newman H. Fisher, Jr. 

SUMMARY 


The shock-expansion method for the calculation of the pressure dis- 
tribution on cylindrical wings in supersonic flow is here extended to 
tapered wings made up of single-curved (i.e., developable) surfaces. The 
method applies In regions of the wing tdiere (a) the con^onent of velocity 
normal to the surface rulings is supersonic and (b) the flow is not influ- 
enced by the presence of the root or tip. The calculation is carried out 
by regarding the single-curved surface as the limit of an inscribed poly- 
hedron TThose edges coincide with rulings of the actual surface. The 
changes of flow across an element consisting of an edge and subsequent 
face of the polyhedron are then found from elementary considerations of 
infinitesimal plane waves and simple geometry. The result, in the l i m it 
of the curved surface, is a pair of simultaneous, nonlinear ordinary dif- 
ferential equations for the components of Mach number normal and tangen- 
tial to the surface rulings. These equations are readily integrated by 
standard numerical methods in any given case. Calculations are carried 
through in the present report for a biconvex triangular wing of aspect 
ratio 4 at two values of the free-stream Mach nimiber. 

As usual, the small-disturbance assumptions can be used to provide 
considerable sinq)lification. In the linear case, in particula r , a closed 
expression is readily obtained for the pressure distribution on a wing of 
biconvex section. For the calcxilated example, the relationship bet^reen 
the linear and the shock-expansion results is similar to that observed 
in two-dimensional flow. In the hypersonic case it is fOTind that the 
pressTore distribution can be calculated by disregarding the taper and 
treating the streamwise section from a simple two-dimensional point of 
view. Examination of two-dimensional shock-expansion results for the 
calculated exaaq)le sh ovtb, in fact, that a strip approximation of this 
kinii provides good accuracy even at moderate supersonic Mach numbers. 


INTRODUCTION 


The present report is concerned with a type of invlscid supersonic 
flovr that occurs on tapered wings made up of single-curved surfaces. By 
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"single-curved” it is meant that the svirfaces are capable of being devel- 
oped onto a plane. Such developable surfaces fall into three classes: 
cylinders, cones, and convolutes, the last being the relatively uncommon 
surface generated by a straight line moving tangent to a nonplanar curve 
(see, e.g., refs. 1 and 2). The supersonic flow over an infinite cylinder, 
yawed or unyawed, can be treated by the two -dimensional shock-expansion 
method, first proposed by Epstein (ref. 3) and since discussed in numer- 
ous books and papers (see, e.g., rrfs. k- and 5). The extension of the 
shock-expansion method to portions of wings made up of cones - or, if 
desired, convolutes - is the subject of the present paper. 

As in the usual two-dimensional case, it is required in the present 
work that any shock wave associated with the leading edge of the wing be 
attached and that the flow Just downstream of the wave be supersonic nor- 
mal to the edge. If these conditions are met, the method then applies in 
those regions of the wing Tdxere (a) the component of velocity normal to 
the surface rulings is supersonic and (b) the flow is unaffected by the 
presence of tips or Junctures. Such regions, thou^ seemingly rather 
special, are of interest for two reasons. Eijrst, at the hi^er super- 
sonic Mach numbers - say 3 and above - they may constitute a majority of 
the wing surface. Second, as will be seen below, the flow in these 
regions is of a simple type which is amenable to a relatively precise 
analysis . In particular, the flow quantities are (with a minor qualifi- 
cation in the cane of the convolute) constant along the surface rulings. 

The reasoning on "vdiich the foregoing statement is based is as fol- 
Icnfs: Consider that the surface of the wing is replaced by an inscribed 

polyhedron whose edges coincide with rulings of the actual surface.^ In 
view of the single-curved nature of the surface, this can obviously be 
done. If the conditions of the preceding paragraph are satisfied, the 
flow at the leading edge of the wing will be characterized by an oblique 
shock wave or a yawed Prandtl-Meyer expansion. In either case the flow 
q.uantities will be constant on the face of the polyhedron adjoining the 
leading edge. At subsequent edges of the polyhedron -vdiere the component ’ 
of velocity normal to the edge is supersonic, an expansion fan will orig- 
inate. If one neglects the reflected effects that arise when these 
expansion fans interact with a leading-edge shock wave,^ the fans them- 
selves will be of the yawed Prandtl-lfeyer iype. To this approximation, 
therefore, the flo\7’ quantities will be constant on each face of the poly- 
hedron. Now let the number of faces of the polyhedron increase indefi- 
nitely. In the limit the polyhedron will become the single-curved surface, 
and the faces of the polyhedron will become the surface rulings. It fol- 
lows that, if reflected effects are neglected, the flow quantities must 

^It will be assumed throu^out that the surface is convex to the flow 
and free of sharp comers, thoxi^ these restrictions are not essential. 

^This approximation is also made in the two-dimensioneOL case. The 
resulting errors are discussed by Eggers, Syvertson, and Eratis in refer- 
ence 6. They conclude that in two-dimensional flow the errors are gmaiT 
except for c ondi tions near shock detachment. There is no reason to expect 
that the situation will be substantially different in the present case. 
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be constant on each ruling, (ih the case of the cone, uiiere the shock 
wave and reflected effects are themselves conical, the qualification in 
this statement can be dropped, and the conclusion taken as exact. It 
can, in fact, be arrived at directly from the familiar argument concern- 
ing the lack of a characteristic length in the boundary-value problem in 
the conical case. For the convolute an argument like the foregoing is 
necessary, and the qualification must be retained. In the calculations 
that follow, as distinct from the argument itself, the reflected disturb- 
ances -vTill always be ignored, and a quantitative approximation 'VTill remain 
in either case . ) 

The method of calculation is arrived at by simple extension of the 
foregoing ideas. To this end, we concentrate attention on an element of 
the polyhedron consisting of one comer and the adjoining downstream face 
and regard the angular deflection at the comer and the angular breadth 
of the face as infinitesimal. Since reflected effects are neglected, 
the expansion occurring at the comer Can now be regarded as taking place 
throu^ an infinitesimal plane wave (i.e., Mach wave) passing throu^ the 
comer . The acconq)anying changes in flow are .readily calculated from the 
fact that the vectorial change in velocity from one side of the ^?ave to 
the other must be normal to the wave. Across the face of the elanent - 
that is, from one comer to the next - the only changes that need be con- 
sidered are those associated with the change in relative orientation 
betvreen the velocity (^diich is constant on the face) and the comers 
themselves. These changes, Tdiich affect the upstream conditions for the 
wave at the next comer, are found from simple considerations of geometry. 
FoUovring this procedure, one obtains, in the limit of the curved surface, 
a pair of simultaneous, nonlinear ordinary differential equations for the 
components of Mach number normal and tangential to the surface rulings. 
These equations are to be integrated from the leading edge rearward. The 
Mach number components at the leading edge, which provide the initial 
conditions for the integration, are found by applying the concepts of 
simple-sweep theory to the equations for an oblique shock wave or Erandtl- 
Meyer expansion. The integration itself is readily accomplished by stand- 
ard nvmierlcal means. 

The first section of the paper is concerned mainly -VTlth the deriva- 
tion of the differential equations and with the evaluation of certain 
derivatives required for their application to conical surfaces. The 
second section discusses the simplifications that occur Tdaen the usual 
small-disturbance approximations are introduced, both in the linear super- 
sonic range and at hypersonic speeds. Of particular interest here is the 
result that at hypersonic speeds the method reduces to a simple strip 
theory in the streamwlse direction. In this section of the paper a closed 
expression is also obtained for the pressure distribution on conical f-rLngs 
of biconvex section in the linear case. The paper concludes with seme 
numerical results for a triangular wing made up of conieeQ. surfaces 'VTith 
apex at the tip. A realization that the pressures on a portion of this 
wing must be constant on rays throu^ the tip was, in fact, the starting 
point of the analysis. 
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Althoii^ the present problem Is a fairly obvious one, reference to 
it in existing literature is difficult to find. The only paper the 
authors have come upon, in fact, is that of Thomson and Sheppard (ref. 7 ) 
fdiich became available after the present analysis vas complete. In this 
reference, the surface of a conical iring is approximated by a polyhedron 
and the pressure distribution on this approximating surface calculated 
by application of simple -sweep theory and the Prandtl-Meyer relations at 
each successive comer. The present plan of incorporating a limiting 
process may be preferable with regard to accuracy and ease of computation 


NOTATICai 
Primary Symbols 


a speed of somd 

c Tdng chord (measured at arbitrary spanwlse station) 

Cp pressure coefficient (see eq.. (30)) 

f(g/c) function defining shape of wing section 

M Mach number vector 

M magnittide of Mach number vector 

p static pressure 

p^ total pressure 

r radial distance along sirrface ruling 

t thietness of wing section 

T absolute temperature 

V magnitude of dlstvirbance velocily, V - Voo 

V veloci-ty vector 


^or the calculation of the entire nonlinear flow field around a 
general conical surface, reference should be made to the work of Maslen 
(ref. 8). The same problem has also been discussed by Perri at the Con- 
ference on Hi^ ^eed Aerodynamics, Polytechnic Institute of Brooklyn, 
January 20-22, 1933 • 
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V 
6V 

x,y,z 
a 
On 
7 • 

8n 

A 

An 

Mn 

V 

g 

5a 

6t 

da 

du 

dr 

du 

<p,u 


magnitude of velocity vector 
increment in V 
Cartesian coordinates 

angle of attack measured in streamwise direction 
angle of attack measured in plane normal to leading edge 
ratio of specific heats (7/5 for air) 

deflection angle measvired in plane normal to leading edge 

local Cartesian coordinates in plane normal to surface ruling 
(see sketch (d)) 

leading -edge angle measured in streamwise direction 

leading-edge angle measured in plane normal to leading edge 

Mach angle corresponding to camponent of velocity normal to sur- 
face ruling 

angle betvreen velocity vector and surface ruling (see sketch (a)) 

chordwise position aft of leading edge at arbitrary spanwise 
station 

infinitesimal angle between two edges of face of polyhedron 
(see sketch (b)) 

infinitesimal deflection angle at comer of polyhedron measured 
in plane normal to comer (see sketch (b)) 

rate of rotation, at a given ruling, of the orthogonal projection 
of a moving surface ruling onto a fixed tangent plane through 
the given ruling 


rate of rotation, at a given ruling, of a moving tangent plane, 
measured in a plane normal to the ruling 

angles used to specify position of rulings on conical surface 
(see sketch (c)) 


Subscripts 


conditions in free stream 
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* conditions at critical speed 

0 values at origin of coordinate system (see sketch (d)) 

1 values at leading edge 

2 values at trailing edge 

n component of velocity or Mach number normal to surface ruling 

(see sketch (a)) 

t con^jQnent of velocity or Mach number tangential to surface .ruling 

(see sketch (a)) 

C change across comer of polyhedron (see sketch (b)) 

F change across face of polyhedron (see sketch (b)) 

Superscripts 

( )* first derivative of function -with respect to argument 
( )** second derivative of function vlth respect to argument 

SHOCK-EXPAHSIQN METHOD 

Differential Equations for Flow Over Single-Curved Surface 


known (sketch (a)). 


Since the flow along a solid boundary must be everywhere tangent to 

the boundary, the velocity vector Y at any point on a single-curved 
svirface ^d-ll be cooqpletely defined if the magnitudes Vn and Vt of its 
components normal and tangential to the surface ruling at that point are 

In the present case, where V is constant on each 

ruling, Vn end V-fc are given 
by a pair of ordineiry differ- 
ential eq.uatlonB that can be 
derived as follows: 

As explained in the 
introduction, we begin by 
regaiding the single-curved 
svirface as the limit of a 
polyhedron whose edges coin- 
cide with rulings of the 
actual svirface. Let us then 
consider the flow over an 



Sketch (a) 
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infinitesimal element of the 
polyhedron (sketch (h)) con- 
sisting of a comer C and 
the adjoining downstream 
face F. The infinitesimal 
deflection at the comer, 
measured in a plane normal to 
the comer, is denoted hy 
6t. The infinitesimal angle 
between the two edges of the 
face F is represented by 
6a, (The positive directions 
of 6t and 6a are shown in 
the sketch.) The changes in 
Vn and Vt across the complete 
element C plus F will be 
considered as congposed of the 
changes across the comer C 
followed by the changes across 
the face F. 



Under the conditions discussed in the introduction, the changes across 
the corner C will take place throu^ a planar Mach wave originating at 
the corner. Let V, Vn, and Vt here denote the magnitudes of the veloc- 
ity and its normal and tangential cong)onents iq)stream of the comer. 
According to the properties of infinitesimal plane waves, the angle between 
the Mach wave and the upstream face of the comer is then given by 



where a is the speed of sound corresponding to V. The situation dom- 
stream of the comer (see sketch (b)) is completely determined by (a) the 

property that the vectorial change in velocity 6V(j across the comer 
must be normal to the wave and (b) the boundaiy condition that the result- 
ant velocity ^ + downstream of the comer must be parallel to the 

face F. From requirement (a) it follows at once that the change in Vt 
across the comer is zero, that is, 

6Vt^ = 0 (2) 

The corresponding change in Vn can be found by considering the triangle 

formed, by the vectors Vn, 8Vn«(s8VQ}, and Vn + BVn^. From this triangle, 
■whose geometry is dictated by requirements (a) and (b), it foUovrs that 



■8 


NACA TN 3^99 


Vn + 5V] 




Vn 


^ 1^) 

(I - Mn - St) 


With the aid of equation (l), this can be written, to the first order in 
infinitesimal quantities. 


8Vnc = 


Vn 


6t 


(3) 


- 1 


This is, of course, identical to the corresponding equation for two- 
dimensional Prandtl-Meyer flow (cf. ref. 9) except for the appearance of 
Vn iihere V irould ordinarily be. 

The changes in Vn and V-t across the face P are found from the 
fact that the velocity vector is constant in magnitude and direction on 
the face. If v is the angle between this vector and the surface ruling 
(see sketch (a)), we have in general 


Vn = V sin V 
V-t = V cos V 

Differentiation of these expressions gives, for constant V, 


dVn = V cos V dv = V-tdv 
dV^ = -V sin V dv = -Vndv 




Since the direction of the velocity vector does not vary on F, the change 
in V from one edge to the other is due solely to the difference in ori- 
entation of the edges. Application of equations (4) to the calculation 
of the changes in Vn and V-^ across F thus gives to the first order 
(and with due regard for the positive direction of 6 ct) 


(5) 


8Vnj. S -Vt8a 
8V.fcp = (Vn + 8Vn^)8(T = 


( 6 ) 
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The changes in flow across the complete element C plus F can now 
be found by combining equations (2) and ( 6 ), and (3) and ( 5 ). In this 
manner we obtain, in the limit as the infinitesimal quantities tend to 
zero, the following ordinary differential equations for the calculation 
of the flovr over a single-ciarved surface: 


dVn 



dr - V-^do 


(7a) 


dVt = Vndo 


(7b) 


For conq)utatiQnal work it is convenient to rewrite the above equa- 
tions in dimensionless form in terms of the normal and tangentieil Mach 
numbers and Mfc, where % - Vn/a and Mfc = Vt/a.'* To do this the 
defining expressions for and are first differentiated to obtain 


dVn 

da 

( 8 a) 

dJ^=_ 

-ISn-r 

II 


( 8 b) 


The speed of sound a is known to be related to the Mach number by the 
adiabatic energy equation, which can be written (cf. ref. 10 , eq. ( 32 b)) 


where a* is the constant critical speed. Differentiation of this rela- 
tion gives 


da 7-1 

® ^ 1 + 2L^ (Mji 2 + Mfc2) 


■^The analysis could be simplified at this point by making the speeds 
dimensionless throu^ division by the constant critical speed a* Instead 
of by the variable speed of sound a. This leads, however, to sll^tly 
increased labor in the eventual numerical work. 
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The desired differentieQ. equations are now found by substituting this 
expression and equations ( 7 ) into equations (8) and solving for dl^ 
and This gives finally 





dr - Mtda 


(9a) 


dMt 


7-1 

2 - 1 


dr + %dff 


(9b) 


The differential equations ( 9 ) provide a means for calculating the 
Mach number (and hence the pressure) distribution on a single-curved sur- 
face under the conditions outlined in the introduction. To integrate the 
equations for a given surface, the differentials dr and do must be 
revTTitten in the form dr = (dT/dA)dA and do = (dq/dA)dA, where A is 
any independent variable by which the position of the surface rulings can 
be conveniently defined. In the ll^t of the limiting process, the deriv- 
atives dr/dA and da/dA can be given simple geometricsil interpretations 
in terms of the rulings and tangent planes of the curved surface. The 
derivative dr/dA at a given aruling can be identified as the rate of 
rotation of a moving tangent plane, evaluated at the ruling in question 
and measured in a plane normal to the niling. The derivative do/dA is 
the corresponding rate of rotation of the orthogonal projection of a mov- 
ing s-urfane ruling onto the (fixed) tangent plane through the given ruling. 
Expressions for these derivatives for a general conical s\orface are given 
(for a suitable choice of A) in the next section of the paper. The actual 
integration can be carried out by any one of a number of standard, step- 
by-step nvauerical procedures, depending on the accuracy required (see, 
e.g., refs. 11 and 12). Obviously, the integration breaks down when 

% ^ 1 . 


It is important to note that the following sign conventions are 
implicit in the derivation of equations (9): (a) is positive iflhen 

^ is in the direction of positive ds (see sketcjh (b)); (b) M^ is 

positive ^dien ^ is directed av?ay from the intersection of the surface 
rulings. These conventions must be carefully observed in carrying out 
the numerical integration. 

Equations ( 9 ) are, of course, applicable to the swept cylinder as 
well as the cone and the convolute. In this case da is identically 
zero, and Mq becomes independent of M;^ (®®® (9a)). This is as it 

should be according to the concepts of simple-sweep theory (ref. 13 ). 
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For sucli a protlem, however, the methods of the present paper are si: 5 )er- 
flnous, since a solution is readily found from tabular results for two- 
dimensional flow (see, e.g., ref. l4). 


Application to Concial Surfaces 


In the following section equations ( 9 ) will he specialized to a con- 
ical surface, since this is the case most commonly encountered in prac- 
tice, For such a surface, the position of the rulings - and hence the 
shape of the cone - can he specified hy an equation of the form 


9 = 9(“) 


where 9 and' u are as indicated in 
sketch (c). Here x, y, and z are 
Cartesian coordinates fixed in the 
■wing, "With origin 0 at the apex of 
the cone,® The x,z plane is taken 
as the chord plane of the wing. The 
velocity of the undisturhed stream is 
assumed parallel to the x,y plane. 
With w as the independen-fc variable 
(replacing A of the preceding sec- 
tion), equations ( 9 ) can he -written 


are as indicated in 



Sketch (c) 


dl^ = 




- 1 


dr ds 

Mfc — dio 

du d(i) 


(11a) 


dMt = 


— — + — du 


®This apex may be located at either the tip or fonTardmost point of 
the \T±Dig. In -the case of a -wing -with a truncated tip, the apex may be 
virtual rather than actual. 
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The task now is the purely geometrical one of expressing dr/du and do/du 
in terms of the knoi«i function cp(u). 

To do this let us consider a representative surface ruling at Wq, cpo* 
At a point o located a radial distance ro along this ruling (see 
sketch (d)), we establish a rectangular coordinate system in a 


y 



plane normal to the ruling . The ^ direction is taken parallel to the 
x,z plane mth ^ positive in the direction of increasing u; T) is 
positive in the direction of increasing 9 . The trace of the conical 
surface in the plane is then given by an equation of the form 


Ti = n(U 


Expressions for t and a are now needed for differentiation. Recall- 
ing the geometrical Interpretation of the derivative of t following 
equations (9), we can wite (with due regard for the positive direction 
of 6t as shown in sketch (b)) 



T = tan 


( 12 ) 
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>diere the subscript o denotes evaluation at the point o. To i^rite an 
expression for a, we consider an arbitrary point a on the trace of the 
conical surface. Let the point b be the orthogonal projection of a 
onto the trace (in the plane) of the tangent plane to the conical 

surface at the ruling oO. In the li^t of the previous interpretation 
of the derivative of a, we can then write 


a = 


tan 


ob 


An expression for ob can be obtained by considering the planar figure 
obac, where c is the orthogonal projection of a onto the ^ axis. 
Since oc = 5 and ac = tj, we have 



The expression for a can thus be written 



Expressions (12 ) and (l3) inust now be differentiated with respect to 
u, and the derivatives evaluated at the point o. Since tj and ^ will 
be expressed subsequently as functions of a, the resulting expressions 
iTill be put in terms of derivatives with respect to this quantity. We 
thus obtain 



o 
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' (l4b) 


Expressions for ^ and as functions of w can be found by return 
ing to sketch (d) . The points e, f , and g are first established as the 
orthogonal projections of o, c, and a onto the x,z plane. The point 
d is then established on the line ag by drawing cd parallel to fg. 
From the geometey of the resulting figure, it is apparent that ef = oc = 
and fg = cd = ac sin tpo = t] sin (po* We can therefore write, from consid- 
eration of the quadrilateral efgO, 


tan(ci) - Uo) = = 




eO - fg rocos cpQ - ■q sin cpo 


(15) 


and 


sin(u - Uq) = 


ef 




(16) 


An expression for gO is found by observation of the triangle gaO. 
Since ad = ac cos cp^ = q cos <Po dg = cf = oe = roSin cp^, it follows 
that 


tan (p 


dg + ad roSin <Po + 'll cos (Pq 
gO gD 


(IT) 


and hence by combination with equation (l6) that 


sin(u - Uo) 

tan cp jCgSin 9o *Po 


(3B) 


Equations (I 5 ) and (iB) constitute a pair of simultaneous equations for 
^ and q. Solution of these equations gives 


sin(u - Uo) 

cos (PqCos(u - Uq) + sin cp^tan 9 


(19a) 
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T] = ro 


tan 9 - tan (pocos((o - cjo) 
cos(w - (1)0) + tan 9 


(19b) 


Since the relation (lO) will be known in any given case^ these equations 
give ^ and tj . as functions of «. 

Equations (I9) must now be differentiated to find the quantities 
required in equations . This operation gives finally 


= rocos 9o 


(4) o ' 

(5) o ' 'Po 

(1;^)^ = ^o(sixi q>oCos cPo + 9o*’ ) 


where the primes denote differentiation of 9 with respect to its argu- 
ment. Substitution of these expressions into equations (l 4 ) gives, after 
the subscripts are dropped, the following relations for the required geo- 
metrical derivatives; 



dr 

dh) 


cos “9 


( 20 a) 




dcr 

du 


(20b) 
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With the aid of equations (20), the differential equations (ll) can be 
integrated numerically for any conical t/ing for •which cp is given as a 
function of u.® 

On most conical "wings encountered in practice, the apex of the cone 
t- mi be located at the 'wing tip. In this case the contour of the surface 
is usually defined by specifying the shape of the strearavri.se airfoil sec- 
tion - that is, by an equation of the form 


- = -<-) 
c c \c/ 


■where c is the chord of the wing at some arbitrary spanwise station z, 

I is the chordvrf.se distance aft from the leading edge at that station, 
t/c is the thickness ratio of the airfoil section, and f is some given 
function of |/c. In an application of -this "type, it is advantageous to 
replace du in equations (ll) by [dw/d(|/c)]d(|/c) and integrate directly 
vrfth respect to |/c. The necessary expression for dco/d(|/c) can be found 
by reference to sketch (e) . If toi and ug are the 'values of w corre- 

, sponding to the leading and trailing edges 

. of the ■wing, it foUovre from the sketch 


that 


^ _/ 1 1 \ 

^ “ ^^tan *" tan uy 

^ ^ y 1 1 ^ 

\tan Ui tan Ug/ 


(22a) 


(22b) 


Sketch (e) 


ELimination of z between these equations 
and solution for tan u gives 


tan 0 ) = 


tan 

I / tan ci)i> 
” c y " tan Way 


®Dr. Stephen H. Maslen of the WAGA Lewis Laboratory has pointed out 
to the au'fchors that equations (U) and (20) can also be derived from the 
general differential equations for conical flow as given in reference 8. 
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from whiclx it follows that 


d(g/c) 


sin^O) 
tan 0)i 



tan oji 
tan 


(lizJ 


(24) 


To find cp as a function of S/c, reference is made to sketch (c). This 
gives 


tan (p 


^ 

z/sin 


u 


\diich becomes, with the aid of equations (21 ) and (22b), 

t sin (0 / tan toA /|\ , . 

tancp=-T (1-T (25) 

^ c tan (jix \ tan tjgy \cy ' 

Equations (20), (23), (24), and (25) provide the relations necessary to 
integrate equations (11 ) as functions of |/c. 


Initial Conditions 


The initial conditions for the integration of equations (9) or (U) 
are the values of the Mach nvnnber con5>onents at the leading edge. These 
values, denoted by Mn^ and 

are found by appl^ng 
the concepts of sin^jle- sweep 
theory (ref. I3) to the 
locally cylindrical flow at 
the edge. 

The situation is as 
shovm in sketch (f).^ The 
given quantities are the 
magnitude M^o of the free- 
stream Mach nuniber, the angle 
of attack a, the angular 
location ui of the leading 
edge, and the leading- edge 
angle 1\ measured parallel 
to the x,y plane. To 



Sketch(f ) 


^ T’To understand this sketch it should be recalled that the vector 
Mqo Is parallel to the x,y plane. 
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apply sin5>le-En7eep theory, one must know the magnitudes and of 

the coorponents of free- stream Mach number normal and tangential to the 
leading edge together with the deflection angle measured in a plane 

normal to the edge. As can be seen from the figure, the Mach number 
con5>onents are given by 


Mfc^ = MooCOs a cos(jt - Ui) 
= -Modcos a cos Ui 


(26a) 


and 


Mn^= n/Moo^ - Mt„2 

= Moo^/l - COS^ COS^Ul 


(26b) 


The deflection angle (taken as positive when the deflection causes a com- 
pression) is given by 


&n = ^n “ o>n 


(27) 


where and on are the leading-edge angle and the angle of attack, both 
measured in a plane normal to the edge. These angles are related to the 
given angles by the following equations (readily derived on the basis of 
the sketch) : 


an = tan~ ^^^ ^ (28a) 

\sin Wi/ 

■Kn = (aSb) 


With Mn^ and 8n known, the value of Mn^^ is easily found from the 
well-known results for an oblique shock wave (8n > O) or Prandtl-Meyer 
eaqDansion (Bn < O) . (See, e.g., refs. 5, 10, or I5.) As required by 
sin5)le-sweep theory, the calculations are carried out in the same manner 
as for two-dimensional flow (see ref. 15)^ except that "normal*^ quantities 
are used in place of the usual resultant values. The value of Mtj^ is 
found from the fact, given by simple- sweep theory, that Vt at the lead- 
ing edge is the same as in the free stream, that is. 
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Vtx = Vt^ 


Division of this equation by the speed of sound gives 


Vtoo Vteo Sco 
ax ax &0O 

or finally 


Mti 


Mtoo 

(ax/a») 



■ (29) 


The value of is given by equation (26a) . The value of or 

Tx/Too is found from two-dimensional theory as before (see ref. I 5 ) . 

As a general remark, it may be noted that the angle of attack enters 
the entire scheme of calculation only throu^ its influence on the initial 
conditions. The differential equations ( 9 ) and (U), and the derivatives 
that appear therein, are all independent of a. 


PressTire Coefficient 


The value of the local pressure coefficient 

'30) 

is easily found once the distributions of Mn, and Mt are known. To find 
p/pj^ - and hence Cp - we distinguish two cases as follows: 

bn < 0.- In this case there is an esjpansion at the leading edge and 
the flow is isentrqpic. The total pressure is therefore everywhere the 
same, and we can write 


P. ^ (P/Pt) 

Poo (P/Pt)oo 


(31) 
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The values of (p/p^)pg and (p/p^) are found from the usual conrpressihle- 
flo^r tables (e. g. , ref. lO) as functions of Moo aiJ^d M, vhere 
M = ^/Mn^ + Mfc2. 

> 0.- In this case^ because of the increase in entropy throu^ 

the leading-edge shock, it cannot be said that p^ = p^ . We can, however, 
write that p^ = p^^ and hence that “ 


Poo (P/Pt)x W/ 


(32) 


The values of (p/p^)^^ and (p/p-j;) can be obtained from the conapresBible- 
flow tables as functions of the corre^onding Mach number. The value of 
(Pi/Poo) is found, as before, from an oblique-shock calculation based on 
the flo\7 normal to the leading edge. 


aiALL-DISTDRBMCE APPROXIMATION 


It is of interest to examine the form the analysis takes when the 
small-disturbance approximations are Introduced. To this end we define 
a dist\rrbance speed v according to the relation 


V = Voo + V (33) 

and assume that terms of hi^er than the first order in v/Voo can be 
neglected. Consistent with this approach, we also assume that the wing 
surface y = y(x,z) lies sufficiently close to the x,z plane that only 
first-order terms in y and its derivatives need be retdlned. 


GenereQ. SmaJl-Di sttirbance Equation 


To incorporate the foregoing assumptions into the analysis, we return 
to the differential equations ( 7 ) and rewrite them in terms of the result- 
ant velocity V. This is done by means of the relation 


= Vn^ + Vt^ 


which becomes after differentiation 
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V dV = Vn dVn + Vt dVt 

By substitution of the expressions for dVjj^ and dV-t from equations (7), 
the foUo^Ting differential equation for V ib obtained: 



The speed of sound a is related to V by the adiabatic energy equation, 
which can be witten 



(35) 


Now let us sin^jlify equation ( 3 ^) in the li^t of the previous assump 
tions. To begin, it can be shown that to a first order in small quanti- 
ties 


dT 


1 

sin ( 1 ) 



Here u again denotes the azimuthal angle of the surface ruling (cf . 
sketch (c)), thou^ the considerations are not now restricted to a cone. 
To the same accuracy, the normal velocity Vn can be written 


^n “ VgaSin 0 ) + (terms of first order) 

where the form of the first-order terms is immaterial for present purposes 
To the first order in v/Va,, equation (35) also sin 5 )lifieB to 

a^ = a«,^|^l - (7 - l)Moo® ^ 

If these three expressions are substituted into equation (3lj-) and the 
assunption is made that MooSin oj » 1 (i.e., the component of free-stream 
Mach number normal to the rulings is not near l) , one then obtains to the 
first order 
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sin oi 


2 2 
Moo sin (i) 


- 1 


1 - (7 - l)Moo^:^ 
Vqo 



(36) 


This eq-uation is general in the sense that it contains both the linearized- 
supersonic and hypersonic approximations. These two approximations "vrill 
now be discussed separately. 


Linear Approximation 


If Moo is not too large, we can write that (7 
equation (36) slnplifies Immediately to 


l)Moo^ ^ « 1, and 

Voo 




sin (i) 


JKci^bId.^ - 1 



(37) 


This resiILt is linear in that the distributions of v/Voo over two single- 
cui^ed surfaces of identical plan form - say yi(x,z) and ya(x,z) - will 
upon being added equal the distribution of v/Voo over a third stirface 
ys(x,z) where ys = yx + ys* is easily verified that equation (37) 
conforms to the similarity rules of linearized sraper sonic wing theory. 

Any velocity distribution calcTilated by means of this equation must, in 
fact, be identical to that which would be obtained by solution of the par- 
tial differential equation of linear theory under the same boundary con- 
ditions. This follows from the fact that the only phenomena i^ecifically 
neg^Lected in the present development - that is, the disturbances reflected 
from a leading-edge shock wave - disappear in the conventional linear 
development as well. Equation (37) will be integrated for a specific 
type of ■VJing (i.e., for a specific relation between 3 y/Sx and u) follow- 
ing a discussion of the hypersonic approximation. 


Hypersonic Approximation 


When Moo is very large, it can no longer be said that 
(7 - 1 )Moo^ ^ « 1 . In this case, however, the first term under the radi' 
cal in equation (36) will predominate, and the equation can be reduced to 
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Ji-W- ^ 


The dependence on u thus disappears from, the differential equation. 
Equation (38) is, in fact, the same equation one obtains if the hyper- 
sonic small-disturhance approximation is applied to the differential equa- 
tion for two-dimensional Prandtl-Meyer flow.s The same situation can he 
shown to exist with regard to the calculation of the initial conditions 
and the pressure coefficient - that is, the equations are independent of 
0) and are the same as those obtained from the two-dimensional small- 
disturbance theory of hypersonic flow. It thus appears that, for thin 
wings at hi^ Mach numbers, the pressure distribution can be found by 
disregarding the taper and treating the wing section from a simple two- 
dimensional point of view. This result is in agreement with the ideas 
expressed by Eggers in reference l6,- It is a consequence, of course, of 
the fact that when the surface Mach mmiber is everywhere large, as will 
be the case on a thin wing at a large free-stream Mach number, the region 
of dependence of a given point on the wing surface is of negligible extent 
in the spanwlse direction. 

Since (0 does not appear in equation ( 38 )> the integration can be 
carried out at once. This will not be done here, since the results can 
be shown to agree with the equations for thin two-dimensional airfoils 
at hypersonic speeds as worked out completely by Llnnell (ref. I7) . These 
equations are arrived at by Linnell by making suitable approximations in 
the usual two-dimensional shock-expansion relations. According to present 
considerations, these approximate two-dimensional results may be applied 
directly to the streamwlse sections of sufficiently thin tapered wings. 

For engineering applications to thicker wings, one mi^t go a step farther 
and apply the complete shock-e:q)ansion method in the same manner. This 
latter possibility >rill be examined in the subsequent example. 



Application of Linear Approximation 
to Conical Biconvex Wing 


The linear approximation (37) will now be applied to a conical wing 
with apex of the cone located at the wing tip. For such a wing the sur- 
face is given by equation (2l) as 


®It can be shown that equation (38) conforms to the hypersonic simi- 
larity rules. 
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This can he differentiated to obtain 


^ ^ ^(y/c) ^ 3(y/e) ^ t yr\ 

Sx 3(x/c) S(|/c) c \c/ 


and 




t 

c 


f • /!) 5! 

\oy au/a(i/c) 


Where the primes denote differentiation of f with respect to its argu- 
ment, Substitution of the last relation into equation (37) and insertion 
of dti)/d(|/c) from eqmtion {2k) then gives for the governing differential 
equation 


a(V\ , . M a.. 

tan (dg 

Since in linearized wing theory the pressure coefficient is given by 
Cp = -2v/Voo, one can also write 


= W ^ _ t^ : 

tan U 2 


If is the pressure coefficient on the upper, surface at the leading 

edge, it follows that the value of Cp at any other point on the same 
surface is given by 




u 


.IT 


(l/c) 


sin 0 ) i^Meo^sin^ - 1 


du 


( 39 ) 


tan U 2 
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IDhe value of Cp^ can be found by applying slnrple-sweep concepts to the 
known results for linearized two-dimensional flow. This gives 




v/Woo^sin^i - 1 ^ 


(J*o) 


On a wing of biconvex section, f**(|/c) is constant. (This is exactly 
true if the section is made tip of parabolic arcs and true to the accuracy 
of small-disturbance theory if it is made up of circular arcs.) Equa- 
tion (39) can then be written 


Cp - Cp^ + 



f**tan (1)1 
^ tan a>i 
tan u)s 


L 


u 


du 


sin (ji - 1 


The integral in this equation can be put into a standard form by means of 
the substitution p = sin^. In this manner one obtains finally the fol- 
lowing equation for the pressure coefficient on a conically tapered wing 
of biconvex section: 


Cp = Cp^ + 


(t/c)f ’ *tan Wx 
^ tan o)x 
tan a>2 


I (Moo^ + l)sin^ - 2 

|sin 

i (Moo^ - l)sin^ 


, -1 (Moo^ + l)sin^u>x - 2 
sin ^ 

(Moo^ - l)sin^fc>i 


m 


The value of Cp is given by eqiiation (40); u is related to the more 
common variable ^S/c by equation (23). Equation (4l) could be obtained, 
presimably, by the conventional methods of linearized conical-flow theory, 
though the authors are not aware that this has been done. 
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KDMERICAL EXAMPLE 


Calculations have heen carried out on the basis of the foregoing 
methods for a sweptback triangular "wing of a^ect ratio 4 (sketch (g)). 

The surface of the wing is tapered 
linearly from the root to the tip - 
that is, each half of the wing is a 
cone with apex at the wing tip. With 
the origin of coordinates at the ri^t- 
hand tip, it follo\?s that Ui = 45° 
and 0)2 = 90 °* The wing section, which 
is symmetrical about the chord line, 
is 5 percent thick with maxl.Tmnn thick- 
ness at midch’ord. It is composed of 
a circular arc from the leading edge 
to the midohord and a parabolic arc 
fraa the midchord to the trailing 
edge. The trailing edge is blunt 

t/c*.05 with a thickness one-half the maximum 

thickness. Q?he function f(|/c) is 
therefore (cf. eq. (2l)) 



Sketch (g) 




for 


0 < l/c < I 


and 


<l) » H - (IJ 


for 1/2 ^ l/c ^1. For this profile f**(|/c) is seen to have a discon- 
tinuity at midchord. 

The calculated values of Cp on the tipper surface of the wing are 
shown in figure 1. Results are given for angles of attack of 0°, ±3°> 
and ±6° at free-stream Mach numbers of 3*36 and 2.k6. In each case, curves 
are shown as ccanputed by the shock-expansion method of this report 
(eqs. ( 11 )), by the corresponding linear approximation (eq. (4l)), and by 
the ordinary (i.e., two-dimensional) shock-expansion method for the stream- 
tfise section. The ntmerical integration of equations (U) was carried out 



NACA TN 3k99 


27 


in the present case by means of the simple "trapezoidal formula" (see, 
e.g., ref. 11 , pp. 2h-2^ or ref. 12 , pp. 236-238). Values given by this 
method ■vri.th ten equal increments of |/c were an5>ly accurate con^jared 
■vrLth results obtained by more refined methods of integration. The curves 
from the two-dimensional shock-expansion method are included in the figure 
for the reasons mentioned in the discussion of the hypersonic approxima- 
tion. Sketches indicating the approximate regions of applicability of the 
results are shown for each Mach number. These regions, as approximated 
on the basis of the free-stream Mach angle, constitute 69 and 56 percent 
of the plan- form area at the respective Mach numbers of 3.36 and 2,k6. 

More accurate determination of the regions as functions of Moo and a 
could, of course, be made on the basis of the results obtained in the 
pressure calculations themselves (cf . ref. j ) . The value of a at which 
Mix becomes sonic just behind the leading-edge wave is -15.2° at Moo = 3*36 
and -6.9 at Moo = 2.1t6. 

The relationship observed in the present case between the shock- 
expansion resvilts (solid curve) and the corresponding linear approximation 
(dashed curve) is much the same as that known for two-dimensional flow 
(cf. ref. 4 , fig. E,3i). As in the latter case, the shock-expansion val- 
ues of Cp are, with minor exceptions, more positive than the linear 
values by an amount which increases as the absolute value of Cp 
increases. As a consequence, the normal force on the wing at a given |a| 
(as indicated by the vertical distance between the Cp curves at equal 
±a) is concentrated more toward the leading edge in the shock- expansion 
results. This effect would, of course, be less ptonounced for wings of 
smaller thickness ratio. 

Judged on the basis of the present shock- expansion findings, the 
curves obtained by the two-dimensional method are a good approximation at 
the positive angles of attack shown in figure 1 . As the angle becomes 
negative, however, the approximation becomes progressively less accurate. 
This is due primarily to the error involved in the calculation of the dis- 
continuous flow changes at the swept leading edge, an error which becomes 
larger as the absolute value of the deflection at the leading edge 
increases. On the present wing, this deflection is zero on the upper sur- 
face at a S 6°. The two-dimensional approximation would be expected to 
deteriorate again at positive angles greater than this value. This behav- 
ior is, of course, dependent on the amount of sweep at the leading edge. 

On a reversed triangular wing, for exanple, where the leading edge is 
uns^rept, the two-dimensional method would give accurate initial conditions 
at all angles, and the over-all situation would undoubtedly be different. 

The logical argument for the use of the two-dimensional method in the 
present case rests, as has been seen, on the assumption that the free- 
stream Mach number is very large. The method would therefore be expected 
to diminish in accuracy as the Mach number is reduced. The results of 
figure 1 shovr this effect. It may be noted, however, that even at a mod- 
erate Mach nTimber such as 3 . 36 , the two-dimensional method provides remark- 
ably good accuracy within the angle range studied here. This agrees -vrith 
previous results (e.g., ref. I8) which show that other deductions based 
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on the hypersonic small-dlstnrhance assun 5 )tions (as, e.g., the similarity 
rules) can he applied with good accuracy at fairly low supersonic Mach 
numbers. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., April 8, 1955 
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